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INTRODUCTION 


The equations for the plane motion of viscous fluids of constant 
volume are, after elimination of the pressure and introduction of the 
stream function ^ hy which the velocity components 


V = V 

X y Bx 


are expressed, reduced to the one equation 


^ ^LMM.= aAA\lf (I) 

at dx ay " ay ax 


therein 0 indicates the ratio between viscosity coefficient and spe- 
cific mass n, eind A signifies the Laplace operator. 

This equation is satisfied by all potential motions 


Alt' = 0 


however, this fact is of little significance since viscous fluids adhere 
to solid walls and, from well-known considerations of function theory, 
there cannot exist a potential motion which wo\ild do so. Otherwise, 
properly speaicing, only Poiseuille's laminar motion is known as exact 
solution of equation (l) and that solution does not even show the sig- 
nificance of the quadratic terms because they identically disappear 
there . 


*"Spiralformige Bewegungen z^er Fliissigkeiten." Jahresber. d. 
deutschen Math. Ver. 25, ISYJ , pp. 34-60. 
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Under these circumstances it seems perhaps usefvil to know a few 
more exact solutions of equation (l) for which the quadratic terms do 
not disappear; such solutions. will be indicated below according to two 
methods . 

In both cases, one deals with motions in spiral-shaped streamlines 
(which are observed frequently) . 

Third, we shall, in addition, investigate the nei^borhood solu- 
tions to pure radial flow. 


FIRST PART 


We raise the question: 

Are there solutions of equation (l) which are not potential motions 
for which, however, the stream paths eure the same as for a potential 
motion whereas the velocity distribution is to be different? 

We shal l be able to indicate such solutions, in fact all of them: 
the streamlines 6ire logarithmic spirals (including concentric circles 
and pure radieil flow); for the velocity distribution, one arrives at an 
ordinary differential equation which for pure radial flow leads to 
elliptic functions. In the discussion, the influence of the quadratic 
terms becomes manifest in a considerable difference between inflow and 
outflow (see paragraphs 7^ 9)- 

We require, therefore, solutions ilr of equation (l) for which 


i|r = f (q>) 


and A9 = 0, but not Ailf = 0. The latter condition excludes 


f«'((p) = 0 


We limit ourselves to steady motions ^ = 0. 

at 

1. The calculation becomes clearer if first the auxiliary problem 
has been solved: 
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Transformation of equation ( I) into isometric coordinates, that is 
such curvilineeu: coordinates <P,X that 

cp + iX = w(x + iy) = w(z) 

Let us thus assume 




^ _ ^X 

Sx by by bx 


If one denotes by there results first, with the abbre 




viation 


dv 

dz 


= Q 


Ai = Q A’f 

With the double integral extended over an arbitrary region, one has 



|A |1 . |1W dy = /T d A d* 
OX oy by bxj Jo 


// 

MQA') 

^ _ 

^(Q A') 


d<P dX 



^x ’ 

bx 

dtp_ 


If 

^(QA>) 

^ _ 

d(Q A') 

bi 



acp 

bx ’ 

dx 

dtp 

\dx dy 

by bx) 


dx dy 


Since, however. 
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there follows 


^ A ^ 
dx by 


b A ^ 
bx 


UY _ U uy _ q2 



M _ 

B A' 


. A In ft St . 

S In ft St^ 

[V5<p 

^x 

^X 

bp) 

^ V Sep Sx 

Sx bp)_ 


However, 


^ In ft _ p 

bv 


bp 


,2 

d w 

In ^ = 2R In ^ = 2R 
dz 


dz 


dw 


m 


and 


^ In ft '- 2 

dx 


Sx . 


R In 


dw 

dz 


-2 A J In ^ = 


bfP 


dz 


-2J In ^ = 


dw 


dz 


-2J 


d^w 

_dzL 

^dw\^ 

<dz/ 


are valid. If one puts the analytic function of z 



= a + hi 


one obtains 


^ A ^ 
bx by 


b A M ^ 
by bx 


, 4 
dw 

fb A* 

M. _ 

S A* 

^Ua- 

faM 


dz 

V ^<P 

Sx 

Sx 

sepy 

V Sx 

Sqv 
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Finally, there results 


A At = Q A' (Q A't) = Q^A' A' t + Q A'Q A't + ^ + 

\ dcp dtp 

= Q^l^i + ^ii|r a - 


a A*t 

hx W 


In Q = In 



is a harmonic ftmction, thiis 


A' In Q = 0 


hence , 


Aia = ( A ln ._af + lA ..af = a2 + b2 

<i \ ) \ J 


Thus one obtains as the result of the conversion of equation (l) to 
isometric coordinates T, X for steady motion 

^ A' t ^ ^ A' -l M + A*t/a M + b = a I A' A't + 

dX bX dT \ ^X ^9/ I 


A>(a3 . b3) . 


la - 


bX ) 


therein, a + hi is the analytic function 

(w = 9 + iX, z = X + iy) 


,2 

d w 

^ 2 
dz 


(g) 

and A' denotes the operator 


(II) 
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2. We return to the question on page 2: ijf must be a mere function 

of «P 


t = f(9) 

If derivatives, with respect to <p, ajre denoted by primes, equation (ll) 
becomes 


f ' 'f 'b 


a + f " (a^ + b^) + 2f ' ' 


(III) 


f may depend only on 9 but must not depend on X. 


This is certainly possible if a and b do not depend on X, thus, 
since a + bi is eui analytic function of 9 + Xi, do not depend on 9 
either, if a + bi is, therefore 


a + bi 



that is, constant. We shall see later (paragraph 3 ) that this is the 
only possibility- 

From a and b being constant, there follows 


w In - zn^ + wn 

a + bi \ ® 

thus, after introduction of the polar coordinates 


z - Zq = re 


di 


9 = - — -(a In r + bd) + 9^ 

a^ + 0 
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Thus , the streamlines 9 = const are identiceil with the logarithmic 
spirals 


a In r + bd = const 

a = 0 signifies pure radial flow, b = 0 flow in concentric circles. 
The velocity distribution, however, is given by equation (ill): the 

rEidial conq>onent is 

^ dJ> 2b f, 1 

r dd r bd g2 ^2 r 

the circular coinponent 

= -f • ^ = 2a 1 

■ dr " Sr ^2 + ^2 r 

consequently f ' ^ the magnitude of the velocity. There- 

fore, f must disappear on solid welHs. 

Without restriction of the generality, one may presuppose left- 
hand spirals so that r increases with d, thus a and b have dif- 
ferent signs; since, furthermore, -(9 + iX) is an analytical function 
just as 9 + iX, and equation (ill) is actually invariant with respect 

to a simultaneous sign^ change of 

a 
b 


Therefore, positive velocity 

f ' > 0 

in contrast for 

f • < 0 


9,a,b, one may presuppose 
> 0 
^0 

components signify for 

r dd dr 

outflow, 

inflow. 


^dJranslator's note: The original says "time change," obviously a 

misprint . 
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Since 9 may be replaced by c9, one may in addition impose a con- 
dition on the constants a and b. 

3> We now want to conduct the proof that on the basis of our require 
ments a and b must be constant, that therefore the flows in loga- 
rithmic spirals are the only ones the flow patterns of which correspond 
to a potential motion without themselves being a potential motion. 

If a and b were not constant, the analytical, function a + bi 
would produce a conformal transformation of the 9 + iX-planej by virtue 
of equation (ill) which with the abbreviations 


A 


^ I ! I 

f ' ' 




(f ' ' = 0 is excluded) may also be written 


a^ + b^ - 2A(9)a - 2B(9)b + C(9) = 0 


(III') 


the circles (equation (ill')) would correspond to the straight lines 
9 = const in this transformation. 

These circles wovild therefore have to form an isometric cxirve 
family. 

However, if the family of curves 

g(a,b,9) = 0 (in') 

is to be an isometric one so that A9 = 0, the function g must satisfy 
the equation 


8,,b®b) ‘ ° 

8a^ 8b^/ 

and this equation must either he identically satisfied, or he a conse- 
quence of equation (ill*). 
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One has 


g^ = 2(a- A), Sb = 2(b-B), Ag = 4, = -2A' , = -2B' 

gjp = -2A'a - 2B'h + C, Scp 9 = 'b + C" 

g^2 + gj^2 ^ 4(a. - A)2 + 4(b - B)^ = 4(a2 + B^ - C) 

Thus, equation (IV) is quadratic in a and b; an easy calculation shows 
the result that the quadratic terms are automat icaJJ-y eliminated. There- 
fore, the coefficients of the two terms must be zero whence follow three 
conditions 


0 


A" C - 2AA' - 2BB' _ B' ' 



C - 2AA' - 2BB' 



C" C - 2AA' - 2BB' 



Hence, there follows further that A’, B’, C, must be proportional and 
fiirthermore that 



must be constant. 

Ihe final, condition yields 

C = cLj^B + p A = r^B - 8 


or with 




f^ = off ' + pf " 


and 


f * * * 


= 7f 'f ' ' + 5f " 
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which, integrated, yields 


f" = i.af'2 + ^f' +€ 


and 


f " = 1 7f '2 + 5f < + 


Comparison of the two values for f'' resxilts in 


- of + Sf + 


_ 2 


Of'^ + pf + 6 
7 f' + 6 


which must he identical with the preceding value of f ' ' . The compari- 
son requires 

7=0 a = 0 p = 6^ e=6q 

thus C = A^ = 6^ constant 6ind f' = 5f + q. 

The second condition 


® = 1l_ = const 

A^ + - C 


= const and this together with f'^ = 6f ' + q 

^ 1 2 

contradiction 


f * = const 


however would yield 
wo\ild result in the 


therewith^ the proof has been produced* 

4. We now turn to the determination of the velocity according to 
differential equation (ill) which may be integrated once and assumes, 
after introduction of the quantity proportional to the velocity at unit 
distance 


u = f(cp) 
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the form 


u’ ' + 2au' + ufa^ + _ -L + C = 0 

' ' - 2a 


This equation is identiceil with a damped oscillation which takes place 
under the influence of the potential 


- ^ u3 + i(a^ + b^)u^ + Cu 
6a 2' ' 

We start with the limiting cases: 

1. The streamlines are concentric circles; b = 0. 

Then 

u =-^ + e'^-'^CA + B9) 

and, because of 


<p =- 2 In j. 


u = const + r^(A + In rj 


whereby, the velocity distribution 


a r 

is given. The exact solution of Conette's case is also contained therein 
the three constants here are determined from the two limiting values of 
the velocity and from the fact that in case of a full turn around the 
circular annulus, the pressure must revert to its initial veJ-ue. An 
easy calculation yields Bq^ = 0 and thus 

- ^ 1 ^) 

(More details on the. determination of the pressure are seen in para- 
graph 10 . ) 
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2. The flow is p\irely radial; a = 0. 
The differential equation reads 


u' ’ = b^u - + C = 0 

2a 


and leads to elliptic functions 


u' = \fToI u3 + 3abu^ + const u + const 
where the three e's are only subject to the one condition 


+ eg + e^ = 3ab 


but otherwise are still at disposal. 

Since,, according to the remark on page 8 one relation between a,b 
is still unused, it will be expedient to put 


b = -2 


so that one obtains, according to page 7 

qp = •a 

Then the conditional equation for the e reads 

®1 ^ ®2 ^ ®2 “ 

and one has 
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thus 


u = -2a + 


-j={i - ^o) i 82-83 

\/6a 


where are the three integration constants. For the pressure 

j 

(see paragraph 10) there results the equation 


^ /P ^ 1 = 1 f f ' ' + f ’ 

Sqp\u 2 j J .2 J .2 


= A 1 A f *2 
3tp j.2\2 


+ 2of 


■) 


its uniqueness is a priori ensured, thus does not determine here any of 
the constants. 


Discussion of the Radial Flow 
5. The condition 


®l ®2 ®3 = 


( 1 ) 


requires at least one- e to have a negatively real part, for instance 
then 

^(^3) “ 

with the equality sign being valid only ^en all three e's have the 
same real part. 

Furthermore, since this part is real, there must apply 
(a) for three real e's either 

-«» < u ^ e^ ^ -2a 
or 

eg S u ^ ei 
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(b) for one real e 

-00 ^ u < e 


where, however, this e may be positive. 

Furthermore, two possible types of flow must be distinguished: 

1 . Either there are no solid walls, thus a source or sink in an 
unlimited fluid. Then u must be a periodic function of 9 , with a 
period which is an integral part of 2n. u = -00 is excluded, u = 0 
need not occur. Therefore, this case can occur only for three real e's> 
and 

e£ ^ u ^ e^^ 

must be valid. 

2 . Or there are two solid walls, for instance for ■d =0 and for 
d = (which may also be equal 2rt); then at these walls u must be 

u = 0. 

(a) In case of three real e's there must be, additionally 

®2 = ° > 0 

and either 

(a) ^2 - ® 

or 

(^) 0 ^ u ^ 

(b) In the case of one real e, this e must be positive and 

0 ^ u ^ e 

One remembers, furthermore, that according to page T, paragraph 2 , 
u > 0 signifies outflow, u < 0 signifies inflow; so that one has 
inflow in the case of 2(a) (a), and outflow in the case of 2(a) (p), and 
2 (b) above. For the case 1 , both cases may occur. 
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^At the distance r = ! 


15 
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thus 


a f 

« Jn 




° sjl + yj: 


II . I 


where 


0 e < 1 


Thus the relation (2') between Ujjj,&,n,a reads, due to the significance 


of 






= 2 
n 6o 


or 


6a + 3% - I y 


( 2 ") 


with Tj being a proper fraction. 
Since, furthermore 


Jo 


% dt 


\jlTx^si^ 


Pz X 

^0 r 


dt 




2^2 


= sinh"^^ I 


thus becomes arbitrarily large with increasing %, one has lim e 

%—oa 

thus, 

^ ■ 


= 0 , 


lim TJ = 1 

X=oo 
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From equation (2'') there follows 

which, with u = 1, gives as the minimum value 

/ 

The mean inflow velocity is therefore considerably limited upward, the 
more so, the easier movable the fluid. 

However, this is the only restriction: If Uj^^ and the integer n 

are selected so that 

6a + 3Uj^ > 6a — 

there exists, certainly, a pertaining 8. 

1 12 
For if — 8 increases from zero to the value 6a + 3%i7 x ^ ^ 

2 “2 

2 

lies between zero and 6a + 30 ^ (because for the second value 9C 

2 12 
becomes infinite and, hence, tj = l) so that certainly sometime ^ H ^ 

2 

becomes equal to 6a + 3Uj^ - 6a — ^daich is presupposed to be positive. 

One sees, furthermore, that for a prescribed period number n the 
fluctuation 8 and for a prescribed fluctuation 8 the period number n 
must increase to infinity with the mean velocity u^. 

Second Case: Outflow Between Solid Walls 

7. The cases 2(a)(3), said 2(b) may be summarized thus 



du 


- u)(u^ + 2au + b) 
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e > 6 is the maximum velocity (at the distance r = 1); because of 


2a = -e^ - e^ = 6a + e 


and p 




> 0, 


otherwise, however, arbitrary 


+ 2au + p 


may for prescribed e assiame all values from + 2oai to oo, so that 


di = 


■ -#/; 


du 


\J{e - u)(u^ + 2au + p) 
appears not at all restricted downward, but upward restricted by 


since 


one has 




du 

u)u(u + e + 6a) 






_ 2j( I 3^ 

a,max \/2e(l + c) + 12a 


(M 


where € signifies a positive proper fraction. 

For the outflow, the width of the wall opening appears therefore 
restricted, according to the preceding equation, by the maximum value e 
of the velocity. For small velocity and large viscosity, the majcimum 
lies near «, otherwise, however, lower; with increasing e it drops 
below all limits. 
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If, therefore, an angle opening smaller than n is prescribed, it 
permits an outflow only up to a certain maximum value. If a greater 
outflow quantity is prescribed, the jet will, therefore, actually prob- . 
ably separate from the walls. 

Also, there is, of course, for any prescribed angle a flow 

possible \diere partly inflow partly outflow occurs. 


Third Case: Inflow Between Solid Wa3J.s 

8. There remains the case 2(a) (a) 

e„ ^ u ^ 0 
2 ~ 



where 

2a - + 6^) “ ^ ®2 

B = e e < 0 

13 


otherwise, however, arbitrary. Thus, the angle -3^^ may be male arbi- 
trarily RTtia 1 1 for prescribed e^. On the other hand, however, it may 
also be made arbitrarily large: one takes, for prescribed e^, the 
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negative value e^ sufficiently close to as far as this is not made 

impossible by e^ < 0. The sole relation between the e 

+ eg + e^ = -6a 

however, resiilts with e^ > 0 in 

-eg - e^ > 6a 

If -Cg ^ 3cr, may actually be assumed arbitrarily close to e^. 

If the maximum inflow velocity is larger than any angle 

is possible between the solid walls* 

If, however, -eg. < 3^, say €3 ct, where € is a positive proper 
fraction, only 

-e^ = e^ + (6 - 3e)a = -eg + + 6(l - e)a 

is possible and 



attains its highest value for e^ = 0 
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where and 6 are positive proper fractions. Thus, the maximum of 
is larger than n. ^ 


When the maximum inflow velocity is smaller than 3^, othe angle 
openings of the solid walls also may attain any magnitude up to «. 

Flow in Spirals 

9* Because of the dajiqping 2au' (see paragraph 4, page 11), a 
periodic solution, aside from u = const, is, not possible. 

A free motion in logarithmic spirals is always a potential motion. 
In contrast, there exist other flows on logarithmic spirals between 
solid walls. 

In order to have, for r = const, the variable cp agree with the 
angle d, one may fturthermore prescribe for the constants a,b, in such 
a manner that one obtains 




2b 


2 

a 


+ b 


= 1 



thus 

b = -1 ± yjl. - a^ 


a must be a proper fraction, otherwise it remains arbitrary. 
Equation III, once integrated, (see page 10) then reads 


u' ' + 2au' + ^ + C = 0 

40 


? I 2 2 

where P = -2b = 2 + 2\jl - a < If, but > a . 
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The velocity at unit distance is 



u is, therefore, the velocity at the distance 



If one first omits the dancing, one has exactly the same case one 
had before except that 



is in front of the square root (see page 12). The relation for the e , 

remains the former one. Since 3^ < 4, the angle opening is increased 
by this influence 

• The damping, however, takes effect in the same sense. Nevertheless 
the main result remains correct. 

For outflow the admissible angle opening is restricted by the 

Tnayimum flow velocity in such a manner that it tends toward zero ^en 
this velocity increases beyond all limits. 

If one puts 


u 


ve 


-a9 


the above differential equation becomes 


V 


I I 


+ 




+ £ v^e"®'*’’ + Ce®'*’ = 0 
4a 
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If <p = 0 is assianed to be the location of the maxinrum Vq for v , 
multiplication by 2v' and integration yields 


* {f? . a^)(v2 - ^ v^a-^^dv a 2C f a=^dv = 0 

From the corresponding equation for u 

u'^ ^ J' u' du + - Uq^) + + 2C(u - iiq) = 0 

one can see that for equal Uq the u-curve will be the steeper, 
therefore the, smaller, the larger C . For value close to Uq this is 

immediat^';^ clear from the differential equation, for in case of u' = 0, 
u" vkuT'be the smaller, the larger C, thus |u'| the larger; from 
the preceding equation one may see, however, that for larger C 


u'^ + 4a 


^^0 


du = u'^-4a/ u'du 

J u 


has the hi^er value. Hence, follows directly for the ascending branch 
(u' > 0) that always |u'| has the hi^er value when C is the larger 
value. For if u' would once reach for the initially flatter curve 

/u 

u' 


du would have to 


have for the former the smaller, thus 


J 

^ u 


du the hi^er value 


which for u' > 0 immediately leads to contradiction since, up to then, 
that is between u and Uq, u' had been the s ma ller value. If, how- 
ever, u' < 0, one has in case of a variation of the C by ^ 


u'^.. + 1» aIu'I du = 

“o'" "o’ "Ju 


= 2 AC 
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or, since 


= A|u'|(2|u'| + A|u' I) 


A|u'|(2|u'| -t-A|u'|) ^ 4a 


dn + 2 AC 


If there were at one point a|u'| = 0, then at this point the first 
term wovild, for fixed AC, decrease with decreasing u, that is, go over 
from positive to negative values; the second term also would decrease 

Lq 

a|u'| du supervening with 

decreasing u wovild. he negative. This is irapossihle, however, since 
the sum of both terms is supposed to be constant 2 AC. 

Therewith, it has been generally proved that the angle opening 
decreases with increasing C (for fixed Uq) ; since the maximum possible 
is desired, the minimum admissible value for C may be assumed. 

O 

This veilue is determined from' v' >0 


since the part of the integral 


I 


2c r° Av s 

Jv 




2 -a9 
V e 


dv 


whence, one may see that the minimum admissible value of C is zero or 
negative^. 

The above inequality must be valid for all v*s between Vq and 
zero and for the positive and negative 9 attained.. One may write them 

2C ^ -(3^ - + v)e"^^° - + VqV + v^)e ^ ^ ^ 


According to page 22, 


p2 _ 


> 0 . 
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where ^Q^cpg are certain mean values. One must further note that, for 
V = Vrtf l^r.1 and |*Po I must be zero vdiereas they have max i m u m values 


V V Q ^ |^ 0 | I 2 I 

for V = 0. The severest restriction is due to the absolutely smaller 
value of the right side, thus the minimum admissible C is given by 


U2 2 \ -“”0' 2 -“(” e '%’) 

= .(p -a)voa -^''oa 


cPo*,^2* positive and the inaximum of 9 q and 9^ which occur for 

V = 0. 

Consequently, one has for the maximum possible -6^ 
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and since 



whence follows 



(Compare formula (4), pa^e 19)- 


Hence follows that with increasing v^, thus also with increasing 
Uq, ^1 drop below al 1 limits: a certain width of the spiral 

permits only a limited outflow velocity. 


SECOND PART 


10. Since the only spiral motion, possible without walls, of the 
type used so far, lead to be a potential motion, exact steady and non- 
steady two-dimensional motions in free spirals will be investigated 
according to smother method. 

In polar coordinates the differential equation (l) reads 


^ Ai + l( h Ai ^ . a At = o A Ai 
^ ry ^r hrj 


A 


r 



r2 ^2 


where 
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I 

Obviously this equation permits solutions which are linear in ^ 


\|; = u + 


in order to make the velocity which has the components 


and v„ = - — - cp ^2^ 

^ ^r Sr 


unique and thus enable a free motion, % must be constant. By this 
statement, the differential equation becomes 



with 


^ LM = a Ai^kl 
St ^ Sr 



Here it is necessary also to investigate the pressure lest perhaps 
in a motion about the singular point r = 0 a multivaluedness of the 
pressiire becomes evident. 

Now one may write the equations of motion without elimination of 
the pressure in the following form 



== At d^ + 




or because of the invariance of the last term 
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/ Hence , follows 



At ^ - r ^la A - AU = Jl 


Au 


or 


^ Au 


+ r 


dt dr 


By virtue ^of the differential equation for u the right side is con- 
stant; thus it must be zero to make the pressure in case of a revolxrtion 
about r = 0 revert to its former value, so that one obtains 


2 

r ^ ^ ^ ^ * Au = 0 

dr dr dt ® 


which by introduction of 



assTsnes the form 


d^v _ A + 2^1 dv _ dv 
^^2 \ ajr ^ dt 


(V) 


Steady Motions 

11. The solution independent of t is 


« + 2 


u = c^r 


+ Cg In r + c^ 


when — < -2, otherwise, \rtien — = 
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If one disregards the trivleil case of potential motion, a spiral motion 
the velocity of which disappears at . infinity exists when 


^ + 1 < 0 
a 


that is, X < - a, thus a sufficiently strong inflow taikes place. 
The spirals then have the form 

« + 2 

cp =-iu = C^r*" + Cg In r + 


If ^ + Z < 0, they approach at infinity the logarithmic spirals; near 
a 

the sink, in contrast, they converge considerably less pronouncedly 
toward the sink point, and the vortex velocity is considerably higher 
than in case of potential flow in logarithmic spirals. 


Unsteady Motions 
12. If one uses the formulation 


V = e“^X (r) 

n 

one obtains from equation (V) 



thus, with the abbreviation X = 1 + — 
^ 2a 
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where the J eure the Bessel functions 


m 


J I / _ — r \ = const 

At 


1 + 5 :^ 


a 4 1(1 + X.) 


* • 

fj2\2/ 2'(i + x)(2 + \) 


If X, does not happen to be an integer, ~X 

regarded as independent solutions. 


13. Sindlarly to the case of the heat conduction equation there 
exist also of equation (V) integrals which show for r = 0 and t = 0 
an indeterminate point. 

Since the differential equation (V) remains imchanged if v is 
multiplied by an arbitrary factor, r by a similar factor, t by its 
square, there must exist solutions of the form 


= “ r°tPw(z) 


After substitution, one obtains for w the differential equation 


, + y./g + 1 - x\ + y/gg. + i\ = 0 

[ ^ j [ kz^ 7 


(VI) 


When does this equation permit a solution of the form 


w = e»»^ 


A slnq)le calcilation yields 


^ = -1 
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6ind then either 


or 


<x = 0 3 = \- l 


a = 2\ 3 = -X - 1 


One thiis has two sln^>le integrals of the required type 


= t^-le 


JL 

4q 


r 

t 


and 


V 


r2X|.-l-Xg 


1^0 t 


for X = 0, both etre transformed into the known integral of the heat 
conduction equation. 

Let us continue the discussion of the differential equation (Vl). 

The singular point z = 0 is a determinate point. The determi nin g 
equation reads 


+ p(a - X) + ~ = 0 

4 


and has the roots 



2 
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so that generally there exist developments of the form 


x-i, 




z + c z + . . . 
2 


6Uld 


a 


Wg = z ^^1 + c^'z + Cg’z^ + . . 


thus 


2)uP-^‘"t 

= r"^H 


2 _> 
'1 ‘'s 


1 + Ct + 


ihct) 


and 


V = t 
2 




1 + C/ + C ■ + 


with the power series continuously converging since z = 0 is the only 
singular point of the differential equation. 

If one assumes p = 0, that is, if one desires solution of equa- 
tion (V) of the form 



an in^.egration by definite integreils is possible. 

The differential equation (Vl) reads after Introduction of the 
roots P^^Pg 


^2 ^ + 2 Ji(l - p - p + z 
,2 dz\ 1 2 

dz 


) + ^^1^2 = ° 


(VI') 
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The connection with Gauss' equation for the hypergeometric fxmction can 
be easily recognized. If one maJces the Euler transformation 


w 



n 

y(s) ds 


with the integral extended over a suitable closed path, one finds for y 
a differential equation which may be satisfied for 


and for 


Therefore 


and 



ty 




by 



'' ^ M - if 


are integrals of equation (VI'). The integrals are extended best over 
a path which leads from R(s) = +» around the points s = 0 and s = z 
back to R(s) = +<». 



Since 
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is einalytically regular in the neighhorhood of z = 0 , there is 


and 



w 


2 



Zi 


-P2 -1+p 
s "^ds 


One can show that 



are the general solutions of equation (v) and likewise are represented 
by definite integrals in closed form. I shall perhaps refer back to 
this and to the connection with the representation and the development 
in terms of Bessel f\mctions elsewhere. 


THIRD PART 

Nei^borhood Solutions to Radial Flow 


l 4 . We shall first look for steady neighborhood solutions to the 
radial flow (pages 12 and I3) by putting 

^lf = f(tp) + p(cp,r) 


where p is assumed to be a small quantity, the square of which is 
neglected. 

We then obtain for f the former equation 


f (IV) + I I + 2 ^ Q 

a 
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with f = u ajid. integrating once 


u' ' + l*u + — + C = 0, 

a 


For p we obtain 


where 


^ Ap ^ u ^ Ap _ 2 u' 

^ r ^r ll- 39 

r 


^ ^ = a A Ap 


3 3r 




+ 


r2 ^cp2 


- ‘ • Since the differential equation in the steady case remains unchanged 
when r and p each are multiplied by an arbitrary factor, there must 
exist solutions of the form 


p = 1^(9) 


One obtains for w the differential equation 


w^ + w' ' (ZX^ -4X + ^ + — u-^u\ + ^ u'w* + 

\ a a / a 

. U3 + 2\2 + u + ^ u* =0 (VII) 

We are partic ular ly interested in free flows and thus in periodic solu- 
tions in 9. 

As concerns the imiqueness of the pressure (see paragraph 10, 
page 28 ), one obtains by a slaq>le calctilation, the condition 

J' - (a. - 2)^d9 = 0 
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On the other hand, there follows from the above differential equation 
itself, by integration over the interval from 0 to 2it with assiuiq>- 
tion of periodicity 


- 2) r wju - (X - 2)o] 

^0 


dcp = 0 


so that in general the uniqueness of the pressure follows from the peri- 
odicity except for the case when X = 2 . 

For free flow, u itself is a periodic fimction of <P; the period 
is an integral w part of 2rt- However, it is not necessary that w 
have the same period as u; but this period must likewise be an integral 
part of 2jt. 

Since 


u’ ' = -C - 4u - i u^ 
a 

•as well as 

u'2 = -2Cu - 4 u 2 - ^ u3 + D = - u)(u - e^) (u - 

may be rationally ejqpressed by u, it will be useful to introduce u 
instead of <P as independent variable in equation (VIl) . Because of 


w' = — u* 
du 


w' ' = ^ ^ u'^ + ^ u' ' 
du2 du 


= ^ u'3 + 3 


du- 


d w 
du^ 


u'u' ' + — u 
du 


11 2 2 
^IV _ d w ^ g d_w u*2^i I + 2 u' +4 u'u' ' ' + ^ u^ 


du 


du'' 


du 


du 


du 
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because of 



all coefficients of the new equation axe integral and rational in u; 
indicating the degree, one writes them 


R. ^ + R. ^ + R. ^ + R- ^ + R_w = 0 (VII' ) 

^,3 ^.2 3du 2 

du du"^ du 

with / 

Bg = u'** = - u)2(u - e^f(u - 


From the fom (VIl) one can see that w possesses singularities 
only where they occur for u, thus . certainly not in the real part of 9 
(\^ich is of interest); ecjuation (VII*) shows that^ as a fimction of u, 
w becomes singular only at the branch points e^^Cg^e^- 


Since = 6u*^u* * is divisible by (e^^ - " ^ 2 ){^ " ^ 3 )^ 

the points e^^e^^e^ are determinate points*, and since the degree of 

the coefficients decreases steadily by 1 with the order of the deriva- 
tives, u = 00 also is a determinate point; the differential equa- 
tion (VII*) belongs to the Fuchs class. 

A well-known calculation yields as the four roots of the determining 
equation for the points , e the values 



Althou^, therefore, two root differences here are integral, no loga- 
rithmics appear in the developments: For from the form (VIl) there fol- 

Ipws that at the points ^ for which u becomes = e, where, therefore, 
u and u* . are regular functions of 9, w also must be such a func- 
tion, whereas ln(u - e) does not possess this regular character. 
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Therefore, the solutions of equation (VII') have at every point e the 
form 


u = Pj^(u - e) + yju - ePgCu - e) 

' other singularities do not exist in a finite domain. 

For u = 00 there results the determining equation 


(2p2 + . 3) + I + I ^ 0 


which has the roots 


= l M2=-| 

which are independent of X, and the roots 

^3 


which 8Lre dependent on X. 

/ Conttnuation 

15. Solutions with the real period 2n (this period must he present 
at least in case of free flow) will exist only for certain X. In 
analogy with Hermite's method for Laine's differential equation, one can 
proceed as follows: 

If Wj^,W 2 ,w^,Wj^ are a fundamental system of equation (VIl) , the 
w(<p + 2jt) are expressed homogeneously linearly by the w 





4 

Wy(q) + 2rt) = 




(v = 1,2,3,M 



1^0 
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There certainly exist periodic functions of the second kind, that 
is, there exist solutions w for which 

w(<p + 2fl) = ow(<p) 

This a is a root of the equation of the foiirth degree 


^ - a 

^12 

^3 

^ik 

®21 

^22 - 

®-23 

®-24 

^31 

^32 

a33 - a 

^34 


%2 

%3 

®44 " 


If a periodic solution is to exist, a = 1 must be a root, suad one 
obtains for \ the equation 

D(l,X.) =0 


The characteristic exponents vhich were calculated suggest the 
attenqpts of putting 


w = u + const. 



and w 



^3) 


Eleaientary ceO-Culation yields the following particular solutions: 

1. The trivial possibility w = u for X, = 0 

2. w = u for X = 2, that is 



'If = f(cp) + ar^f'((p) 
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where a must be small and therefore with the same approximation 

t = f (tp + ar^^ 


so that the streamlines are approximately the spirals 

9 = - ar^ 

0 

f ' remains the same elliptic function discussed before in the case of 
radial flow. It is true that this flow now cannot exist as free flow, 
since this is precisely the exceptional case \ = 2 (see page 37) ^ and 
the condition for the uniqueness of the pressure can certainly not be 
satisfied for w = u. 

3- w = u + 30 , when A. = 1 and C = 3^ whence for ^ 

no contradiction results. 


for A. = -1 and e^ = 0. This solution has a 

/ 

u; likewise, w = i ^ ^ X = -1 and e^ = 0. 


5*. w = - e^) or w = - ^2)(^ " ^ 3 ) X = 1 

and ® This solution too has a period twice that of u. 

The large X may he easily calculated approximately from equa- 
tion (VIl). For such large X there is in first approximation 

w^ + 2A,^w' ' + X^w = 0 


4. w = - eg 

period twice that of 


that is, w = (we restrict o\n:selves to the periodic solutions). 


so that ^ is the period, 
approximately integral. 


The large set -apart X-v8J.ues are therefore 


Finally, one case may be ceJ-culated quite elementarily: the case 
when u is constant, the basic flow therefore an all around vuiiformly 
distributed flow. 
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This case is also of significance for the more general one since 

I4. 

according to a well-knovn theorem by Cauchy and Boltzmann the period 
of w in first approximation is obtained if the constant mean value is 
inserted for the periodic u, under the presupposition that the larger 
fluctuation e^ - e„ be sufficiently small. 

J- ^ . . f . 

For constant u there follows from equation (VII), page 36 


^IV ^ I ^ ~ 2i + w^X,^ - + 4X^ + ^ ~ ^ u^ = 0 


thus with the formulation 

w = e^^*^ 


|j,^ - p.2/2X^ - 2X + 4 + Z-=-k - 4x3 + 4x2 + y2 2 , 7 ._X ^ = 0 


This equation has four roots 


= ±X 1^ = (X - 2)2 - 2 - ^ u 

(5 


SO that all integral positive and negative \ are possible (potential 
motions) as well as all \ which are calculated from 



with integral 

For the case u = const, that is: 

16. For the radieG- flow ^ich is uniform all aroxind, the tinsteady 
nei^borhood solutions also can be given. 

^Boltzmann, Ges. Abh., Bd. 1, S. 43- 
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The differential equation now reads (see page 36) 


dt r Sr 


One may integrate it either by means of the formulation 

Xt+ni9 / ^ 

Ap = e w(r) 

(n integral) thus arrives at the differential equation 



which may be solved by Bessel functions , or by means of the foimulatlon 
(conq)are page 32) 


Ap = = e“i'Pr®w(z) 


whereby one obtains for w(z) the differential equation 


z^w' ' + zw 




For z = 0 this equation hsis a determinate point, the determining equa- 
tion has the real roots 


P 


_m_lu + l ^2 ) 1 u^ 
2 ' 4 a ~ 2^ 4 ^2 
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By introduction of the roots and the differential equation 

assumes the form 




+ zw 




p +z)+ppw = 0 
2 / 12 


This is, however, exactly the differential equation (VI') of page 33 so 
that everything said about it there is also valid here. 


Translated by Maxy L. Ifahler 
National Advisory Committee 
for Aeronautics 
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